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Abstract 

Linear  stability  of  the  incompressible  flow  along  a  streamwise  comer  is  studied  by  solving 
the  two-dimensional  eigenvalue  problem  governed  by  partial  differential  equations.  It  is  found  that 
this  fully  three-dimensional  flow  is  subject  to  inviscid  instability  due  to  the  inflectional  nature  of  the 
streamwise  velocity  profile.  The  higher  growth  rates  for  the  inviscid  instability  mode,  which  is  sym¬ 
metric  about  the  comer  bisector,  as  compared  to  the  viscous  Tollmien-Schlichting  instability  opera¬ 
tive  away  from  the  comer  is  consistent  with  the  experimental  findings  that  the  comer  flow  transitions 
to  turbulence  earlier  than  the  two-dimensional  Blasius  flow  away  from  the  comer. 
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1.  INTRODUCTION 


The  last  three  decades  saw  significant  progress  towards  better  understanding  of  the  stability 
of  flow  over  two-dimensional  and  swept  wing  flows.  These  advancements  have  greatly  improved 
our  prediction  and  control  capabilities  of  the  laminar-turbulent  transition  process.  For  further  prog¬ 
ress  it  is  important  to  enhance  our  understanding  of  the  effects  of  geometric  complications  such  as 
wing-body  junction,  finite  wing  span,  and  surface  roughness  elements,  which  play  an  important  role 
in  the  overall  transition  process.  Here  we  consider  the  inviscid  stability  of  the  flow  along  a  stream- 
wise  comer,  which  can  be  considered  as  a  model  for  the  flow  near  a  wing-body  junction.  The  insta¬ 
bility  mechanisms  induced  by  the  streamwise  comer  will  also  help  in  assessing  the  effect  of  side 
walls  on  transitional  and  turbulent  wind  tunnel  experiments.  Understanding  of  the  instability  will 
also  help  us  devise  techniques  for  controlling  transition  in  the  comer  flow. 

The  viscous  flow  along  a  comer  that  is  formed  by  the  intersection  of  two  semi-infinite  per¬ 
pendicular  plates  (figure  1)  is  three  dimensional  close  to  the  comer  due  to  the  strong  interaction  of 
the  boundary  layers  on  the  two  perpendicular  walls.  The  basic  laminar  flow  permits  a  similarity  solu¬ 
tion  under  boundary  layer  assumption.  The  governing  self-similar  boundary  layer  equations  ap¬ 
propriate  for  the  comer  region,  which  blends  with  the  two-dimensional  Blasius  boundary  layer  and 
the  outer  potential  flow  away  from  the  comerline,  was  obtained  by  Rubin  (1966).  Numerical  solu¬ 
tions  to  these  governing  equations  (Rubin  &  Grossman  1971  and  Ghia  1975)  exhibit  a  secondary 
cross-stream  flow  which  is  directed  towards  the  comer  along  the  two  walls  and  directed  away  from 
the  comer  along  the  comer  bisector.  The  resulting  streamwise  velocity  profile  along  the  comer  bisec¬ 
tor  is  inflectional  in  nature.  Therefore  the  three-dimensional  boundary  layer  near  a  streamwise  cor¬ 
ner  is  susceptible  to  inviscid  instability,  while  the  two-dimensional  Blasius  counterpart  is  only  sub¬ 
jected  to  milder  viscous  instability. 

Experiments  on  comer  layer  by  Nomura  (1962),  Barclay  (1973),  El-Gamal  &  Barclay 
(1978)  and  Zamir  &  Young  (1970, 1979)  have  yielded  contrasting  results  for  the  laminar  self-similar 
velocity  profile.  The  differences  among  these  experimental  results  highlight  the  exceptional  sensi¬ 
tivity  of  the  laminar  comer  layer  solution  to  differences  in  the  shape  of  the  leading  edge  and  stream- 
wise  pressure  gradients  (see  Zamir  1981).  This  sensitivity  to  measurements  can  be  attributed  to  the 
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early  instability  of  the  comer  layer  compared  to  the  Blasius  boundary  layer.  Based  on  his  favourable 
pressure  gradient  experiments,  Zamir  (1981)  observes  that  the  zero  pressure  gradient  comer  layer 
becomes  unstable  at  a  Reynolds  number  of  around  104,  while  the  critical  Reynolds  number  for  a  zero 
pressure  gradient  flat  plate  boundary  layer  is  an  order  of  magnitude  higher,  around  105. 

Lakin  &  Hussaini  (1984)  considered  the  stability  of  the  comer  flow  sufficiently  away  from 
the  comerline  in  the  blending  region,  where  the  streamwise  and  wall  normal  velocities  are  given  by 
the  Blasius  solution  with  a  superimposed  secondary  spanwise  velocity  induced  by  the  comer.  Solu¬ 
tions  to  the  stability  equations  were  obtained  with  a  critical  layer  analysis.  Recently,  Dhanak  (1993) 
studied  the  stability  of  the  blending  boundary  layer  numerically,  employing  the  same  governing 
equations  as  those  of  Lakin  &  hussaini  but  emphasized  the  importance  of  enforcing  appropriate  sym¬ 
metry  boundary  conditions  along  the  comer  bisector,  instead  of  the  usual  asymptotic  boundary 
conditions  at  infinity  for  the  one-dimensional  problem.  However,  his  results  close  to  the  comer  re¬ 
gion  are  quantitatively  questionable  since  the  one-dimensional  stability  analysis  employed  in  his 
analysis  ignores  any  spanwise  variation  of  the  mean  flow.  Where  as  the  actual  flow  is  three  dimen¬ 
sional  close  to  the  comer  and  the  approach  towards  the  Blasius  boundary  layer  is  only  algebraic  (Pal 
&  Rubin  1971). 

Here  we  will  retain  the  strong  dependence  of  the  comer  boundary  layer  along  the  two  wall 
normal  directions  and  consider  the  stability  of  this  flow  with  a  locally  parallel  assumption  along  the 
streamwise  direction.  The  resulting  two-dimensional  stability  analysis  results  in  an  eigenvalue  prob¬ 
lems  which  poses  far  greater  computational  challenges  than  the  one-dimensional  counterpart.  To 
simplify  the  analysis  and  computations,  we  will  restrict  attention  to  an  inviscid  analysis  through  an 
extended  Rayleigh  equation  (a  partial  differential  equation).  The  inviscid  analysis  should  be  ade¬ 
quate  to  capture  the  qualitative  features  of  the  dominant  instability  mechanism  arising  from  the  in¬ 
flectional  nature  of  the  streamwise  velocity  component. 

2.  Mathematical  Formulation 

2.1  Mean  Flow 

The  viscous  flow  along  the  comer  that  is  formed  by  the  intersection  of  two  semi-infinite  per¬ 
pendicular  flat  plates  (Figure  1)  can  be  simplified  with  the  boundary  layer  theory.  Unlike  flat  plate 
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and  infinite  wing  geometries,  the  comer  flow  is  three-dimensional  (all  three  velocity  components 
exist  and  they  are  function  of  all  three  coordinates).  Sufficiently  far  way  from  both  the  flat  plates 
(region  I)  the  flow  can  be  modelled  as  potential  flow.  Close  to  the  plates  but  far  away  from  tne  comer 
line  (regions  II  and  III)  the  mean  flow  is  nearly  two-dimensional  and  depends  only  on  x  and  y  coordi¬ 
nates  in  region  II  and  on x  and  z  coordinates  in  region  ID.  These  two-dimensional  blending  boundary 
layers  are  primarily  Blasius  boundary  layers,  but  with  a  superposed  transverse  flow.  In  the  region 
close  to  the  comer  line  (region  IV)  the  coupling  that  is  created  by  the  mutual  interaction  results  in 
a  strongly  three  dimensional  boundary  layer.  The  flow  in  this  comer  region,  termed  as  the  “comer 
layer”  is  the  subject  matter  of  this  paper. 


The  governing  three-dimensional  boundary  layer  equations  appropriate  within  the  comer 
layer  can  be  written  in  a  self-similar  form  as  follows  (Rubin  1966): 

-9#-f#  +  #  +  #  =  o 


drj  *3£  '  drj  '  a£ 

-  rj  _  3£m  +  d^u 

vudTj  a?  Bn1  ap 

-  uv  -  V  u—  -  £  u—  +  M 

u  vudr,  Qud^vdn  +  Wd^  37  + 

-  -  &  + 

1  Btj  3£  drj  as  d£  drj1  as2 


(1.1) 

(1.2) 

(1.3) 

(1.4) 


where  rj  and  £  are  the  nondimensional  boundary  layer  coordinates  along  the  wall  normal  directions 
y  and  z  given  by 

’ '  7f?Re'/2  ?  =  7frRe‘/2  (2) 


where  the  Reynolds  number  is  defined  as  Re  =  .  The  dependent  variables  u,  v,  vv  and  p  are 

the  nondimensional  velocity  components  and  pressure  and  are  related  to  their  dimensional  values 
u,  v  and  w  and  p  through  the  following  relations: 


0  =  ^Rel/2  ,  *  =  ^Rel/2 


(3) 


Here  U  is  the  free  stream  velocity,  x*  is  the  dimensional  distance  from  the  leading  edge,  g  and  v  are 
the  density  and  kinematic  viscosity  of  the  fluid.  The  dependence  on  the  streamwise  direction,  x,  has 
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been  eliminated  in  the  above  equation  (1)  by  assuming  a  self  similar  solution  for  the  comer  layer 
along  the  streamwise  direction. 

The  elliptic  nature  of  the  governing  equation  (1)  requires  boundary  conditions  for  u,  v,  vi* 
on  the  four  bounding  lines;  >7  =  0,  £  =  0  ,  ij  =  -*  «  ,  £  =  £max  oo .  The  proper 

boundary  conditions  on  the  two  walls  are  no-slip  and  no-penetration.  The  asymptotic  boundary 
conditions  appropriate  in  the  limit  of  £max  -*■  » ,  and  rj  <  £,  should  blend  with  the  blending 
boundary  layer  (region  II).  By  symmetry  the  same  boundary  condition  applies  in  the  limit  of 
£max  00  •  These  asymptotic  boundary  conditions  can  be  expressed  as  the  following  expansion  in 
inverse  powers  of  distance  from  the  comerline  (Pal  &  Rubin  1971) 


♦  u'(,) + ;f> + “f  > 

Smax  S  max  9  max 

(4.1) 

+ ;f> + ;f> 

’max  ’max 

(4.2) 

X,  x  *  /  \  n  W2(7)  .  W3O7) 

w(t],£max)  ~  Wj(>7)  +  .  +  -2 

9max  Smax 

(4.3) 

V  x  _  a  ,  pxb 7)  P2(V)  ,  Hv) 

Pw>Cmax)  Po\V)  ^  c*2  03 

^max  Smax  £max 

(4.4) 

The  zeroth  order  boundary  conditions,  u0,  v0  and  p0 ,  are  nothing  but  the  Blasius  boundary  layer 
solutions.  The  first  order  streamwise  and  wall  normal  velocities  can  be  shown  to  be  zero  (Pal  &  Ru¬ 
bin  1971)  and  the  only  first  order  effect  of  the  comer  is  to  induce  the  secondary  flow,  h»j  ,  towards 
the  comer  along  the  bottom  wall.  The  higher  order  terms  in  the  above  equation  have  been  obtained 
by  Pal  &  Rubin  (1971)  by  requiring  a  simultaneous  matching  of  the  comer  layer  with  the  outer  po¬ 
tential  flow  as  >7,  £  -*■  to  .  Thus  the  above  equation  (4)  provides  asymptotically  accurate  higher  or¬ 
der  boundary  condition  which  can  be  applied  at  the  outer  boundaries  of  a  computational  domain  that 
has  been  truncated  to  a  large  but  finite  rj  =  rjmi%  and  £  =  £max.  The  above  equation  displays  the 
algebraic  decay  towards  the  Blasius  boundary  layer  as  £  -+  » . 

The  governing  equation  ( 1 )  can  be  solved  with  wall  boundary  conditions  and  asymptotic  out¬ 
er  boundary  conditions  (equation  4)  to  obtain  the  laminar  comer  layer.  This  problem  of  obtaining 
the  comer  mean  flow  has  been  addressed  by  many  authors  (Carrrier  1947,  Pearson  1957,  Rubin  & 
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Grossman  1971,  Desai  &  Mangier  1974  and  Ghia  1975)  with  varying  degrees  of  approximations 
applied  to  the  governing  equations  and  the  boundary  conditions.  Here  we  have  employed  a  spectral 
ADI  technique  in  order  to  solve  the  comer  layer  equations.  The  spectral  discretization  will  provide 
exponential  accuracy  needed  for  accurate  stability  calculations.  Results  obtained  from  this  tech¬ 
nique  compare  favourably  with  those  of  Rubin  &  Grossman  (1971)  and  Ghia  (1975).  Figure  2  shows 
contours  of  streamwise  velocity  and  a  vector  plot  of  the  cross-stream  secondary  velocity.  Figure  3 
shows  the  velocity  profiles  plotted  along  the  comer  bisector.  An  inflection  point  in  the  streamwise 
velocity  profile  at  rj  =  5  ***  2 . 4  is  evident. 

3.  Stability  Analysis 


The  stability  of  the  laminar  base  flow  to  small  perturbations  can  be  investigated  through  the 
standard  linear  stability  analysis.  By  making  a  quasi-parallel  flow  assumption  for  the  base  flow 
along  the  streamwise  direction,  velocity  and  pressure  perturbations  of  the  following  form  (normal 
mode  anastaz)  can  be  superimposed  on  to  the  mean  flow 

up  =  wp(J7,?)exp[t(a§  -  at)]  vp  =  Vp(i/,£)exp[r(a£  -  cur)] 

wp  =  wp(J7,£)exp[r(a£  -  cur)]  pp  =  pp07,£)exp[r(a£  -  <ur)]  L" ' 

where  £  is  the  non-dimensional  streamwise  coordinate  and  Jlx* Re-1/2  and  U  are  used  as  the 
length  and  velocity  scales.  These  perturbation  quantities  will  be  added  to  the  nondimensionalized 

base  flow  ( um  =  u  ,  vm  =  v(2Re)-1^2  ,  wm  =  vv(2  Re)  “ l/2)  and  the  total  velocity  and  pressure 
when  substituted  into  the  Navier-Stokes  equation  and  the  incompressibility  condition  and  linearized 
results  in  the  following  stability  equations 


A  A  .A  dllp  A  dup  A  dum  ,  A  dUm  I  A  f  A  A  ,,  ,  \ 

umi  a  up  +  +  Wm-^r  +  vp~^  +  wp~qT  +  1  a  pP  ~  RUP  =  1  w  UP  (61) 


m - p  '  dr}  m  p  drj  p  yp  R  p 

A  A  A  A  dVp  A  dVfn  I  A  SVtfj  .  p  I,  A  A 

‘  a  +  V" if  +  “'-if  +  -  fr,  =  <  <»  vp 


A  A  A  d\Vp  A  dWp  A  d\Vm  .  A  dwm  ,  dp p  A  A  -  _ 

uml  a  Wp  +  Vm-foj-  +  +  vP~dp  +  Wp~df  +  Hf~  Kp  =  1  m  Wp 


A  ^  3vp  dwp 

•au”+-4+if=0 


where  operator  L  —  —  a2  +  and  R  =  Jl  Re.  In  the  above  temporal  stability  formula- 

dtp 
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tion,  the  input  parameters  to  the  stability  analysis,  a  and  Re,  are  respectively  the  streamwise  wave- 
number  and  Reynolds  number  and  i o  is  the  resulting  complex  eigenvalue  whose  real  part  represents 
the  disturbance  frequency  and  the  imaginary  part  corresponds  to  the  disturbance  growth  rate.  Hence 
if  <Uimag_pan  >0  comer  flow  is  susceptible  to  unstable  small  amplitude  disturbances. 

3.1.  Bisector  Instability 


Before  solving  the  partial  differential  eigenvalue  problem  (equation  6),  we  will  first  consider 
the  effect  of  the  inflectional  profile  by  studying  the  simpler  problem  of  the  one-dimensional  stability 
of  the  velocity  profile  along  the  comer  bisector.  Instead  of  the  cartesian  coordinates  rj  and  £  a  new 
orthogonal  coordinate  system  s  and  t  will  be  considered,  where  s  is  along  the  comer  bisector.  In  the 
transformed  coordinates  a  locally  parallel  flow  assumption  is  made  by  ignoring  the  variations  in  the 
mean  flow  along  the  streamwise  £  and  tangential  t  directions  and  stability  of  this  one-dimensional 
base  flow  to  disturbances  of  the  following  form  is  considered 


ipis) 

‘  wPis)  ' 


exp[t(a£  +  fit  —  tut)] 


The  results  obtained  from  this  bisector  stability  analysis  will  be  presented  below  to  provide  qualita¬ 
tive  understanding  of  the  effect  of  the  inflectional  nature  of  the  base  flow  on  the  overall  stability. 
Figure  4  shows  temporal  growth  rate  (imaginary  part  of  c y)  of  the  disturbance  plotted  against  stream- 
wise  wavenumber  a  at  three  different  Reynolds  numbers,  Re=5000  (case  I);  Re=4.5  X 104  (case  II) 
and  Re=1.25  X 105  (case  III).  Two-dimensional  disturbances  corresponding  to/M)  are  the  most  am¬ 
plified  and  results  corresponding  to  this  case  only  will  be  presented.  The  outer  edge  of  the  computa¬ 
tional  domain  is  chosen  to  be  smax=25  and  is  discretized  by  85  Chebyshev  Gauss-Lobatto  points  and 
asymptotic  boundary  conditions  are  applied  at  the  outer  edge.  The  results  presented  are  well  con¬ 
verged  and  show  insensitivity  to  the  exact  number  of  grid  points,  location  and  nature  of  the  outer 
boundary  condition.  Also  plotted  in  figure  4  is  growth  rate  vs  (Re)1/f2  fora =0.21  (case  I V).  Acritical 
Reynolds  number  of  Recrjt=435.0  is  obtained  and  the  corresponding  critical  streamwise  wavelength, 
acrit=0.2  and  the  critical  frequency,  (<yCru)rcal_part=OT .  This  result  compared  with  the  critical  Re¬ 
ynolds  number  of  9.1  X 104  for  the  Blasius  profile  shows  that  the  inflection  point  induced  by  the 
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streamwise  comer  has  the  potential  to  decrease  the  critical  Reynolds  number  by  as  much  as  two  or¬ 
ders  of  magnitude.  Owing  to  its  inviscid  nature,  the  growth  rates  for  the  bisector  profile  are  also 
much  larger  than  those  of  the  Tollmien-Schiichting  disturbance  in  the  Blasius  boundary  layer.  Figure 
5  shows  the  frequency  (tUreaLpan)  variation  corresponding  to  the  four  cases  discussed  above.  As  ex¬ 
pected  the  frequency  varies  linearly  with  a  for  large  values  of  the  wavenumber  and  is  almost  inde¬ 
pendent  of  the  Reynolds  number. 

3. 2.  The  Two-Dimensional  Inviscid  Eigenvalue  Problem 


With  the  above  encouraging  results  we  will  consider  the  two-dimensional  inviscid  instability 
of  the  comer  flow  in  the  limit  of  Re  -*  oo .  In  this  limit  the  viscous  terms  drops  out  of  equation  (6). 
It  should  also  be  noted  that  in  the  limit  of  infinite  Reynolds  number  the  mean  flow  is  purely  stream- 
wise,  since  in  the  boundary  layer  approximation  the  cross-stream  velocities,  vm  and  vvm ,  are  order 

~P=  smaller  than  the  streamwise  velocity.  Equation  (6)  can  therefore  be  simplified  as: 

/Re 


A  A  A  A  dUm  .  A  dtifft  A 

t  <D  up  +  umi  a  up  +  vp-^  +  Wp-^r  =  - 1  a  pp 

(7.1) 

A  A  A  ^P p 

-  «  ®  V,  +  «„ «  a  - -Jji 

(7.2) 

A  A  A  ^P  P 

-  1  (o  wp  +  um  1  a  wp  =  -  -~ 

(7.3) 

a  3vp  dwp 

,au'+st+  a?  =° 

(7.4) 

- - - -  .  - - - — - - - j  . . . . . .  . . . . 

higher  order  equation  for  the  pressure  perturbation  (Hall  &  Horseman  1990): 

*2PD  _  ^dUrn^Pp  ,  |5  ^Pp  _  jBuJPp  =  <o  rA 
dndn  m  dt  dt  a  Lpi 


2  A 
—  n4,  u 


or  um  p.  +  u. 


an2 


18) 


We  will  employ  a  spectral  methodology  (Canuto  etal.  1988)  for  solving  the  above  eigenva¬ 
lue  problem  with  Chebyshev  discretization  along  the  77  and  £  directions.  With  this  spatial  discreti¬ 
zation  equation  (8)  reduces  to  a  generalized  matrix  eigenvalue  problem  of  the  form 

Aa  _  QJ  n  a 

■  ^_p  ~  a=  Pp  (9) 

With  N  grid  points  along  each  of  the  T]  and  £  directions  the  size  of  the  matrix  eigenvalue  problem 
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is  nearly  N 2  x  N2,  therefore  computational  time  and  memory  places  stringent  limitations  on  the  spa¬ 
tial  resolution.  Here  we  have  employed  the  inviscid  approximation  since  the  discretized  matrix  tor 
the  viscous  eigenvalue  problem  (Eqn.  6)  is  sixteen  fold  larger.  The  following  symmetry  conditions 

Symmetric  Mode  Anti-Svmmetric  Mode 


up(y>Z)  -  up(t,Tj) 

Vp(TJ,t)  =  Wp&Tj) 

PpinX)  =  pp(£,i 7) 


Up(rj,Z)  =  -  uP(£,t}) 

vp{yX)  =  ~  7) 

pp(yX)  =  -  Pp(£,y) 


can  be  used  to  reduce  the  size  of  the  matrix  four  fold,  but  the  eigenvalue  problem  still  remains  formi¬ 
dable. 


In  the  inviscid  limit  the  appropriate  boundary  condition  to  be  applied  on  the  solid  boundaries 
for  the  velocity  eigenfunctions  is  no  penetr  .tion.  From  equations  7.2  and  7.3,  the  corresponding 
pressure  boundary  conditions  on  the  solid  boundaries  are  zero  normal  derivatives.  The  computation¬ 
al  domain  will  be  truncated  to  be  a  Finite  square  and  the  outer  boundary  conditions  will  be  imposed 
at  y  -  >7 max  and  £  =  £max.  The  appropriate  boundary  condition  for  pressure  at  these  outer  bound¬ 
aries  is  not  clear.  But  far  away  from  the  comer  (y  =  ymax,  £  =  £max)  the  streamwise  velocity  profile 
is  not  inflectional  and  will  not  support  inviscid  instability.  Therefore,  it  is  reasonable  to  assume  that 
pp  =  0  at  y  =  ymax  and  £  =  £max.  in  any  case,  the  eigenvalue  problem  was  solved  with  both  Di- 


richlet  ip p  =  0)  and  Neumann 


j 


boundary  conditions  and  the  sensitivity  of  the  sta¬ 


bility  results  to  the  placement  of  the  outer  boundary  at  y  =  ymax  and  £  =  £max  and  to  the  number 
of  grid  points  was  also  considered.  Based  on  these  sensitivity  tests  ^max  =  £ma*  =  25with55Che- 
byshev  Gauss-Lobatto  points  along  both  the  y  and  £  directions  was  found  adequate  to  provide  well 
converged  results.  Both  Dirichlet  and  Neumann  boundary  conditions  for  pressure  at  the  outer 
boundaries  yielded  identical  results  to  5  decimal  places. 


Only  results  for  the  symmetric  disturbance  case  will  be  presented  below.  The  anti-symmetric 
disturbances  did  not  provide  any  growing  solution.  This  result  is  to  be  expected  since  the  symmetric 
disturbances  have  their  peak  value  along  the  corner  bisector,  where  the  base  flow  is  inflectional. 
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Whereas,  the  pressure  and  streamwise  velocity  components  of  the  antisymmetric  disturbance  mode 
are  zero  along  the  comer  bisector.  Although,  there  are  no  growing  anti-symmetric  modes  in  the  in- 
viscid  limit,  we  anticipate  growing  anti-symmetric  (Tollmien-Schlichting  like)  viscous  modes  in  a 
full  viscous  stability  problem,  but  their  growth  rate  will  still  be  smaller  than  the  corresponding  sym¬ 
metric  inviscid  mode. 

In  figure  6  the  growth  rate  of  the  most  unstable  inviscid  instability  mode  is  plotted  as  a  func¬ 
tion  of  the  streamwise  wave  number.  The  continuous  curve  corresponds  to  the  case  where  homoge¬ 
neous  Dirichlet  condition  is  used  at  the  outer  boundaries  for  the  pressure  eigenfunction  and  the  sym¬ 
bols  correspond  to  the  case  of  homogeneous  Neumann  pressure  boundary  condition.  Insensitivity 
of  the  results  to  the  exact  nature  of  the  boundary  condition  is  apparent,  indicating  that 
rj max  —  £max  —  25  is  adequately  away  from  the  comer  region.  Also  plotten  in  this  figure  is  tne 
growth  rate  of  the  second  most  unstable  mode,  whose  maximum  growth  rate  is  nearly  two  and  a  half 
times  smaller  than  that  of  the  most  unstable  mode.  The  corresponding  frequency  variations  are 
shown  in  figure  7  as  a  function  of  the  streamwise  wavenumber.  A  nearly  linear  increase  in  frennpncy 
with  wavenumber  can  be  observed.  The  frequency  corresponding  to  the  second  most  unstable  mode 
is  slightly  smaller  than  that  of  the  most  unstable  mode. 

The  inviscid  instability  mode  with  the  largest  growth  rate  corresponds  to  a  streamwise  wave- 
number  of  a= 0.225.  This  result  compares  well  with  results  obtained  from  the  bisector  instability- 
analysis,  where  a  corresponding  to  the  most  amplified  disturbance  is  0.20  at  Re=5000  and  0.21  at 
Re=4.5  X 104  and  further  increases  with  Reynolds  number.  The  maximum  growth  rate  (imaginary 
part  of  i o )  obtained  from  the  two-dimensional  inviscid  instability  analysis  is  0.004  and  is  much 
smaller  than  those  obtained  from  the  one-dimensional  bisector  instability  analysis.  This  is  because 
the  mean  flow  in  the  two-dimensional  analysis  progressively  becomes  less  inflectional  av.ay  from 
the  comer  bisector  and  the  overall  effect  is  to  reduce  the  growth  rate  in  comparison  with  the  one-di¬ 
mensional  analysis.  The  frequency  (real  part  of  a>)  obtained  from  these  two  instability  analysis 
agrees  very  well.  For  example,  the  non-dimensional  frequency  of  the  most  amplified  two-dimen¬ 
sional  inviscid  mode  is  0.1  which  compares  well  with  (GfireaLpart^O'OSQ  and  0.092  for  the  most  am¬ 
plified  bisector  modes  at  Re=  5000  and  4.5  X  104  ,  respectively. 
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Figures  8a  and  8b  show  the  real  and  imaginary  parts  of  the  most  amplified  pressure  eigen¬ 
function  (a  =  0 . 225;  co  =  0  .  108  +  i  0 . 003975).  Peak  values  of  the  real  part  of  the  pressure  ei¬ 
genfunction  occurs  along  the  corner  bisector  near  r]  =  £  =  2.5  where  the  mean  streamwise  velocity 
profile  is  inflectional.  Although  the  peak  values  of  the  imaginary  part  occur  away  from  the  comer 
bisector  closer  to  the  walls,  the  magnitude  of  the  real  part  dominates  the  imaginary  part.  From  the 
pressure  eigenfunctions  the  corresponding  velocity  eigenfunctions  can  be  evaluated  based  on  equa¬ 
tion  (7).  Figures  9a  and  9b  show  contours  of  the  v-eigenfunction  plotted  in  two  different  ranges.  It 
is  clear  from  this  figure  that  the  v-velocity  of  the  disturbance  rapidly  increases  near  the  critical  layer 
where  the  denominator  i(aum  —  co)  that  occurs  in  the  evaluation  of  the  velocity  eigenfunctions  be¬ 
come  nearly  singular.  But  the  pressure  eigenfunction  is  well  behaved  in  this  critical  layer  region  and 
therefore  the  eigenvalue  computations  are  well  resolved.  The  corresponding  w-eigenfunctions  can 
be  obtained  from  figure  9  based  on  the  even  parity  of  the  w  and  v  eigenfunctions  about  the  corner 
bisector. 

4.  Conclusions 

Finally  we  conclude  with  a  few  comments  on  how  the  above  results  are  relevant  in  explaining 
the  rapid  transition  observed  in  zero  pressure  gradient  corner  flow  experiments.  The  two-dimension¬ 
al  inviscid  instability  analysis,  although  confirms  the  possibility  of  an  inviscid  instability  due  to  the 
inflectional  nature  of  the  mean  streamwise  velocity  does  not  provide  any  clue  as  to  the  critical  Re¬ 
ynolds  number  for  this  mechanism  to  be  active.  On  the  other  hand,  the  bisector  instability  analysis, 
although  it  ignores  all  variations  in  the  mean  flow  away  from  the  bisector,  yields  a  critical  Reynolds 
r  umber  of  435.  In  any  case,  these  results  suggest  a  possible  destabilizing  inviscid  mechanism  active 
starting  from  a  point  close  to  the  leading  edge  as  compared  to  the  corresponding  viscous  instability 
in  the  Blasius  boundary  layer.  These  results  are  consistent  with  the  experimental  observations  of  Za¬ 
mir  (1981)  that  the  corner  flow  becomes  transitional  at  Reynolds  numbers  ^  1 04  compared  with  the 
critical  Reynolds  number  of  ^  9  X  104  for  the  Blasius  boundary  layer.  Further  confirmation  of  the 
present  results  with  a  full  two-dimensional  viscous  stability  analysis  (using  equation  (6))  is  currently 
being  pursued. 
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Figure  Captions 


Figure  1:  Streamwise  comer  flow  geometry  and  the  coordinate  system  used. 

Figure  2:  Contour  plot  of  the  streamwise  component  of  the  mean  velocity  field.  There  are  19  contour 
lines  running  from  0.05  near  the  walls  to  0.95  away  from  the  walls.  Also  shown  is  a  vector  plot  of 
the  cross-stream  velocity  components. 

Figure  3:  Velocity  components  along  the  streamwise  direction  and  along  the  comer  bisector  plotted 
against  distance  from  the  comerline  along  the  comer  bisector  (5).  Due  to  boundary  layer  scaling  the 
actual  velocity  component  along  the  comer  bisector  is  order  (Re)1/2  smaller  than  the  streamwise 
component. 

Figure  4:  Growth  rate  obtained  from  the  bisector  instability  analysis.  Case  I:  growth  rate  vs  a  for 
Re=5000;  case  II:  growth  rate  vs  a  for  Re=4.5  X 104;  case  III:  growth  rate  vs  a  for  Re=1.25  X 105; 
case  IV:  growth  rate  vs  Re  fora=0.21. 

Figure  5:  Frequency  obtained  from  the  bisector  instability  analysis.  Case  I:  Frequency  vs  a  for 
Re=5000;  case  II:  Frequency  vsa  for  Re=4.5  X  lO^case  III:  Frequency  vs  a  forRe=l  .25  X 105;  case 
IV:  Frequency  vs  Re  fora=0.21. 

Figure  6:  Growth  rate  vs  a  obtained  from  the  inviscid  instability  analysis  for  the  first  two  most  unstable 
disturbance  modes.  Slid  line  corresponds  to  the  most  unstable  mode  and  the  dashed  line  corresponds 
to  the  second  most  unstable  mode. 

Figure  7 :  Frequency  vs  a  obtained  from  the  inviscid  instability  analysis  for  the  first  two  most  unstable 
disturbance  modes.  Slid  line  corresponds  to  the  most  unstable  mode  and  the  dashed  line  corresponds 
to  the  second  most  unstable  mode. 

Figure  8:  Contours  of  the  real  (a)  and  the  imaginary  (b)  parts  of  the  pressure  eigenfunction  correspond¬ 
ing  to  the  most  amplified  inviscid  disturbance.  For  the  real  part  there  are  21  contour  lines  ranging 
from  -1 .0  (marked  1 )  to  0.0  (marked  L)  and  for  the  imaginary  part  there  are  2 1  contours  ranging  from 
-0.2  (marked  1)  to  0.1  (marked  L). 

Figure  9:  Contours  of  the  absolute  value  of  the  v- velocity  eigenfunction  corresponding  to  the  most  am¬ 
plified  inviscid  disturbance.  For  clarity  the  contours  are  plotted  in  two  different  ranges.  In  figure  9a 
there  are  19  contours  ranging  from  1.0  (marked  1)  to  19.0  (marked  J)  and  in  figure  9b  there  are  20 
contours  ranging  from  0.1  (marked  1)  to  2.0  (marked  K). 
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